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CHANNELS  THAT  COOPERATIVELY  SERVICE 


A  DATA  STREAM  AND  VOICE  MESSAGES,  II: 
DIFFUSION  APPROXIMATIONS* 


J.  P.  Lehoc2ky 
Carneige-Mel Ion  University 
Pittsburgh,  PA 


D.  P.  Gaver 

Naval  Postgraduate  School 
Monterey,  CA 


1.  INTRODUCTION 

In  this  paper,  the  study  of  the  behavior  of  an  element  of 
a  communication  system  carrying  both  data  and  voice  traffic  is 
continued.  In  a  previous  paper,  Gaver  and  Lehoczky  (1979),  a  fluid 
flow  approximation  was  developed  to  predict  the  characteristics  of 
the  data  queue  length  for  such  a  system.  The  fluid  flow  approach 
applied  to  the  case  in  which  the  data  service  rate,  n,  was  large 

compared  with  the  voice  service  rate,  u ;  for  example,  when 

4 

n/u  "  10  .  In  this  paper,  a  different  approach  is  taken:  a  Wiener 
process  approximation  is  developed.  The  accuracy  of  this  approxi¬ 
mation  depends  on  the  'heavy  traffic"  assumption,  that  is  that 
the  overall  traffic  intensity  should  approach  unity  from  below. 

One  need  not  make  any  assumption  concerning  the  individual 
parameters.  It  follows  that  the  diffusion  approximation  comple¬ 
ments  the  fluid  flow  approximation  to  give  a  more  complete  picture 

/ 

of  the  behavior  of  voice-data  communication  systems. 

Research  in  part  sponsored  by  ONR  at  Naval  Postgraduate  School 
N 00 1 4  fl OW  P.O 0 0 6  7 ,  and  in  part  by  NSR  at  Carnegie-Mel  Ion  University, 
ENG79  05526. 
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The  communication  system  to  be  studied  is  essentially 
a  part  of  the  SENET  network  as  described  by  Coviello  and  Vena 
(1975)  or  Barbacci  and  Oakley  (1976)  .  This  network  employs 
time-slotted  frames.  A  certain  portion  of  each  frame  is  allocated 
to  voice  traffic,  while  any  data  traffic  can  use  all  remaining 
capacity  including  any  left  unused  by  voice.  The  voice  traffic 
cannot  use  any  unused  data  capacity  and  operates  as  a  loss  system. 

We  introduce  probabilistic  assumptions,  conventional  in 
many  queueing  studies.  Voice  traffic  arrives  according  to  a 
Poisson (A)  process,  and  each  voice  customer  has  an  independent 
exponential (p )  service  time.  Data  arrivals  are  governed  by 
an  independent  Poisson  (6)  process  and  exhibit  independent 
exponential (n )  service  times.  A  total  of  c  channels  are  re¬ 
served  for  exclusive  use  of  data,  while  v  channels  can  be  used  by 
both  data  and  voice;  however,  voice  pre-empts  data.  The  assump¬ 
tions  imply  that  voice  can  be  modelled  as  an  M/M/v/v  loss  system, 
and  the  well-known  "Erlang  B"  loss  formula  will  give  the  loss  rate. 
We  focus  our  attention  on  the  behavior  of  the  data  queue  and  seek 
to  develop  expressions  for  the  steady-state  distribution  and  the 
mean  queue  length.  Many  authors  have  studied  such  a  system,  in¬ 
cluding  Halfin  and  Segal  (1972),  Halfin  (1972),  Fischer  and  Harris 
(1976),  Bhat  and  Fischer  (1976),  Fischer  (1977),  Chang  (1977), 
and  Gaver  and  Lehoczky  (1979) .  Only  the  latter  paper  gives 
approximations  valid  for  the  extreme  but  realistic  case  in  which 
n/u  is  large,  and  none  give  a  Wiener  approximation . 
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Rather  than  focus  on  this  simple  version  of  the  problem. 


we  wish  to  provide  a  somewhat  more  general  analysis.  The  generali¬ 
zations  arise  when  one  considers  more  complex  voice  service. 

Two  types  of  generalizations  are  rated  here. 

In  general,  voice  traffic  is  unlike  data  traffic.  A  voice 
communication  is  actually  a  series  of  bursts  or  talkspurts 
separated  by  silence.  It  is  reasonable  to  model  a  conversation 
as  a  two-state  Markov  chain  of  alternating  talk  and  silence.  Ex¬ 
ponential  holding  times  are  appropriate  for  talkspurt  and  silence 
lengths,  see  Brady  (1965).  During  the  periods  of  silence,  the 
voice  channel  is  available  for  data  traffic.  This  extra  channel 
capacity  can  serve  to  allow  for  increased  data  utilization,  or  to 
reduce  the  data  queue  length,  or  both.  The  voice  process  can  be 
described  by  a  Markov  chain  for  which  each  state  provides  a 
description  of  the  number  of  voice  customers  requesting  service 
and  the  number  actually  using  the  channel.  For  each  voice  state 
one  can  compute  the  number  of  channels  available  for  data,  and 
thus  the  data  transition  rates. 

A  second  generalization  allows  for  the  addition  of  extra 
voice  channels  by  reducing  the  quality.  One  strategy  available 
to  provide  adequate  voice  and  data  capacity  is  to  assign  voice 
users  to  certain  acceptable  quality  voice  channels,  say  6  or  8 
KBPS.  When  a  specified  number  of  voice  users  are  in  the  system, 
new  voice  users  are  assianed  to  less  acceptable  channels,  say  with 
2  or  4  KBPS.  In  this  way,  one  has  increased  the  number  of  voice 
users  that  can  be  serviced  at  the  expense  of  the  oualitv  of  the  voice 
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transmissions.  This  situation  can  also  be  incorporated  into 
our  structure.  One  again  defines  voice  states:  each  state 
provides  a  description  of  the  voice  channels  in  use.  One  can 
specify  a  generator  Q  for  the  Markov  process  over  the  states 
and  use  it  to  compute  the  channel  availability  for  data  for 
each  state.  To  include  this  and  other  possible  generalizations, 
we  assume  that  the  voice  process  is  a  continuous  time  Markov 
chain  with  finite  state  space  and  that  it  subordinates  the  data 
queue  length  process.  We  will  work  out  formulas  for  the  special 
case  in  which  the  voice  is  an  M/M/v/v  queueing  system. 
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2. 


Diffusion  Approximation  Approach 

We  assume  that  the  voice  process  (V(t),  t  >  0)  can 
be  described  by  a  Markov  chain  with  state  space  {1,2,...,N}. 
We  let  this  Markov  chain  have  generator  Q  =  (q^),  an  N  x  n 
matrix,  and  stationary  distribution  ii  .  An  important  special 
case  occurs  when  voice  operates  as  an  M/M/v/v  loss  system  in 
which  case 
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Each  voice  state  i  gives  a  number  of  channels  available 


for  data,  say  c^.  From  c^  one  can  compute  the  rate,  r^,  at 
which  the  data  queue  increases  or  decreases.  Here  r^  = 6  -  n^. 
We  let 


R  = 


(2.2) 


rN 

For  stability  one  requires  L^=i  <  ^  •  In  the 

special  case  for  which  Q  is  given  by  (2.1),  these  conditions 

<v 

become 


v  Pv 

T  -rr  ( <$  -  (c  +  v  -  i)n)  <  0 
i=0  1* 


where  data  has  exclusive  use  of  c  channels  and  the  voice  has 
v  channels.  This  condition  can  be  written  p  <  1  where 
p  =  (pd  +  pv(l-q) )/(c +  v) ,  Pd  =  5/n,  Pv  =  A/y  ,  and  q  is  the 
blocking  probability 


q 


We  develop  a  Wiener  process  approximation  under  a  heavy- 

rN  4 

traffic  assumption  that  ri7ri  '  This  entials  finding  an 

appropriate  infinitesimal  mean  and  variance.  Locally  the  data 
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queue  length  process  does  not  resemble  a  Wiener  process.  If 
voice  is  in  state  i,  the  change  in  the  data  queue  is  nearly 
deterministic  and  equal  to  r^dt  rather  than  a  normally  dis¬ 
tributed  random  variable.  Nevertheless,  if  we  consider  a  longer 
period  of  time,  say  ft,  t  +  T]  with  T  large,  then  the  change 
in  queue  length  will  be  more  reasonably  assumed  to  be  normally 
distributed.  The  next  section  contains  the  outline  of  a  proof 
that  a  suitably  sealed  version  of  the  process  converges  to 
a  Wiener  process.  This  proof  does  not  require  the  fluid  flow 
approach  of  n/y  large,  so  a  slightly  more  refined  variance 
term  is  derived.  The  Wiener  process  derived  will  have  a  reflect¬ 
ing  barrier  at  0. 

The  voice  process  subordinates  the  data  process;  see 
Feller  [  ],  pp.  for  the  notion  of  subordination.  If  voice 

is  in  state  V(t) ,  then  the  data  queue  is  increasing  at  rate 
r^t)  over  [0,  t  +  dt].  The  expected  change  over  [t,  t  +  T] 
is  given  by 

t+T  t+T  t+TN 

E  J  rV(s)dS  =  J  E(rV(s))ds  =  J  riP(v(s)=i)ds  . 

For  the  case  of  T  large  we  can  assume  V(t)  has  the  stationary 

distribution,  thus  V(s)  has  distribution  tt  for  all  s.  The 

rN 

expected  change  over  [t,  t  +  T]  thus  becomes  T  £^_^r^7K.  Notice 
that  the  data  queue  is  very  large,  so  the  boundary  is  not  en¬ 
countered.  For  a  Wiener  process,  this  expected  change  would  be 
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mT  again  assuming  that  the  boundary  is  not  encountered.  We 
thus  find 


N 


m 


=  l  r.ir.  =  ir  Rl 
“  ^  11  ^ 


(2.3) 


We  follow  a  similar  method  to  determine  a  .  The  variance 
in  the  change  in  the  data  queue  length  over  [t,  t + T]  is  given 
by 

t  +  T 

Var  J  rV(s)dS  • 

We  again  assume  V(s)  has  distribution  tt,  so  this  variance  is 
equal  to 


Var  / 


V(s) 


ds 


-  e(/T  rv(s)ds)2  -  (E  l  rv(s)ds)2  -  rv(s)>2  -  («1)2  T2 


The  last  term  can  be  computed  to  be 


T  T 


T  T 


El[  rV(s)ds)  -  E  l  l  rV(s)rV(t)dsdt  -/  l  E(rV(s|rV(t))dsdt 


T  T  N 


=  2  /  /  l  r  P  (V(s)  =i)  E(r„,  Jv(s)=i)dsdt 

00  i=l  1 


T  T  N 


N 


=  2  /  /  l  r.TT.  I  r.  P(v(t)=j  |v(s)=i 

0  0  i=l  1  1  j=l  3 


=i)dsdt  . 
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We  let  p..(u)  =  P  ( V  (u)  =  j  I V  (0 )  =i)  and  P(u)  =  (p..(u)). 

*ij  ~  lj 

The  Kolmogorov  forward  equations  give  P(u)  =  exp(gu)  where 

2 

exp(M)  =  I  +  M  +  M  /2\  +  •••  for  a  square  matrix  M.  One  can 
reexpress  the  second  moment  in  matrix  form  as 


T  _  T  t 

E  ( /  rv^s)ds)  =  2  /  /  7tR  exp(Q(t-s))  Rl  dsdt 

T  t 

=  2t:r  /  /  exp(Qu)  dudt  Rl  . 

0  0  ~  ~~ 


In  view  of  the  fact  that  Q  is  singular,  one  can  most 

>%/ 

easily  carry  out  the  required  integration  by  introducing  the 
eigenvalue  decomposition  Q  =  $D¥  with 


a  diagonal  matrix  of  eigenvalues. 


$,  ~  &2'  *'*  '  ' 


¥ 


&2' 
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the  associated  right-  and  left-eigenvectors.  Here  $¥  =  ¥$  =  1. 
Since  Q  is  irreducible  and  finite  6.  has  a  negative  real  part 


for  2  <  i  <  N.  it  follows  that 


T  t  T  t 

/  /  exp  (Q  )  dudt  =  $  /  /  exp(Du)dudt 

0  0  ~u  ~  0  0 


0  .  T 

with  ou  =  0  1(0~1(e  1  -1)).  Since  0^  has  a  negative  real  part 
ou  ~  -T/0^  as  T  For  large  T,  the  second  moment  can  be 

written  approximately  as 


thus  the  variance  is  given  by 


-  T2(7t+R1)2. 
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IT-  m  T  - 

=  f-  ££§  (»  -  ®IU  +  ™  J,  • 


For  our  choice  of 


S  1  =  0 ,  thus 


a  It-  l  t  - 

-  tt  rs  m  "  ml)l  =  -  tt  RS  y 


T  - 

where  y  =  (R  -  ml)l.  We  note  that  ir  y  =  0,  so  if  z  =  S  y, 

T 

then  Sz  =  ^  is  a  consistent  set  of  equations  and  tt  z  =  0. 

In  general  one  can  solve  (2.6)  by  direct  calculation 

of  Q  ;  however,  this  is  rarely  possible.  A  second  approach  is 

T 

to  solve  Qz  =  y  subject  to  tt  z  =  0.  For  Q  of  form  (2.1), 
this  can  be  carried  out.  The  equations  become 


Pv(zi  -  z^^)  -  (i-1)  ~  zi-2*  =  ri-l  ~  m  '  1  <  i  £  v 


-V|V  -  Vl>  •  rn  -  m 


(2.7) 
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TT|M 


We  find 


_  ,  v  k  T.  , 

itTrz  =  -  -  z  +  -  l  ir.r.  I  1-1 
y  0  y  k  k  .L 


■  1  p  it  .  i 

1=1  V  1-1 

k  T . 


=  h  l  Vrk'm)  l 

u  k=l  i=l  pv  i=l 


,  v  T.  ,  v  .  v  T.  - 

=  I  l  _ilL_  l  7T,  (r.-m)  =  i  l  1-1 

II  **  OTT.  L  V  V  II 


y  .**,  p  ir  k  k 

1=1  v  l-l  k=J 


i-1 

(-  I  Mr 


y  .  ^  ott.  .  . 

p  1=1  Hv  l-l  k=0 


k'  k 


i  v  Tt  . 

-  V  — — — 

p  1=1  Vi-1 


,  v-i  t: 

-  1  I  — 

n  ^  n  tt 


y  -  A  P  TT  . 
1=0  pv  l 


We  thus  derive 


,  v-1  TT 
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T.  =  l  it.  (r.-m) 
1  j=0  3  3 


For  the  special  case  considered  r^  =  S  -  (c+v-j)n 
=  n (pd  -  (c+v)  +  j).  One  finds 

v  p3 

m  =  l  n(Pd  -  (c+v)+j)  yr|- 
i=0  J  v 

=  n[pd  -  (c+v)  +  p v ( 1— q)  ] 

where  q  =  p^/v!Sv  is  the  blocking  probability.  Thus 
m  =  -n (c+v)  ( 1— p ) .  Plugging  into  (2.8)  we  find 


-m) ) 


(2.8) 
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We  thus  have  given  explicit  formulas  for  the  infinitesimal 

mean  and  variance  of  a  Wiener  process  with  reflecting  barrier  at  0. 

One  can  therefore  use  standard  results  for  Wiener  processes  with 

reflecting  barriers  to  derive  the  stationary  distribution — it 

2 

will  be  exponential  with  parameter  2m/ o  .  The  stationary 

2  2 

mean  queue  length  will  thus  be  a  /2m;  where  a  and  m  are 
given  by  (2.6)  and  (2.8).  For  example  when  v  =  1  and  c  =  0 
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This  is  developed  under  heavy  traffic  conditions,  i 
For  such  conditions 


,  i.e. ,  p  f 


pd  +  rrV  *  1 

V 


pd  -  ttt 


It  follows  that  the  above  expression  is  also 

n  pvpd _ 

y  (1-p)(1+pv)2 

in  agreement  with  Fischer  (1977)  and  Gaver  and  Lehoczky  (1979) . 

The  characterization  of  the  data  queue  as  a  Wiener  process 
with  reflecting  barrier  allows  one  to  discuss  the  dynamics  of 
the  data  queue.  For  example,  suppose  the  data  queue  is  at 
level  x  at  time  t  where  x  is  large.  One  might  wish  to 
study  the  time  that  elapses  until  the  queue  becomes  empty. 

This  is  the  duration  of  a  busy  period  started  at  x,  and  corresponds 
to  a  first-passage  time  for  a  Wiener  process  which  we  denote 
by  t^.  Moreover  one  might  also  be  interested  in  the  area  beneath 
the  sample  path  until  it  reaches  the  boundary  of  Q,  since 
this  area  gives  the  total  time  waited  by  all  customers  involved 
in  the  busy  period.  We  let  this  area  be  A(x)  and  seek  to 
compute  t^  =  E(Tx)  and  a(x)  =  E(A(x)).  Straightforward 
backward  equation  or  martingale  arguments  give 


16 


(2.10) 


tx  »  -x/m 

2  „2 

is  x  .  a 

a(x)  =  '^  +  ^x 

2 

where  m  is  negative, and  m  and  a  are  given  earlier.  The 
distribution  of  the  first  passage  time  can  also  be  easily  cal¬ 
culated;  however,  its  Laplace  transform  is  most  convenient  to 
find  and  is  given  by 

-st 

E(e  x)  =  exp 


+  2s 


(2.11) 
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4 .  A  Refined  Diffusion  Approximation 

In  this  section  we  continue  the  development  of  a  diffusion 
approximation  for  the  data  queue  length  process.  The  analysis  is 
patterned  after  the  important  semigroup  methods  of  Burman  11979] . 
One  develops  a  sequence  of  Markov  processes  and  studies  the 
behavior  of  their  generators.  One  shows  that  the  generators  con¬ 
verge  and  hence  concludes  that  the  associated  semi-groups  con¬ 
verge.  This  entails  the  convergence  of  the  finite  dimensional 
distributions.  The  limiting  finite  dimensional  distributions 
will  be  those  of  a  certain  Brownian  motion.  The  theory  underlying 
this  approach  is  based  on  theorems  of  Trotter  and  Kato,  Kurtz 
and  Burman;  the  reader  is  referred  to  Burman  (1979)  for  technical 
details.  We  will  illustrate  the  approach  in  the  context  of  the 
special  case  v  =  1,  c  =  0.  The  method  applies  to  the  case  of 
general  c  and  v  but  the  details  are  not  given  here. 

Let  ( V ( t) ,  t  2  be  the  voice  process,  a  Markov 

process  with  state  space  {0,1}.  Let  N(t)  be  the  data  system 

size  at  time  t.  We  study  the  sequence  of  Markov  processes 

<{(X  (t),  V  (t)),  t  >  0)>“  ,  where 
n  n  —  n=l 

Xn(t)  =  N(nt)//n  and  Vn(t)  =  V(nt)  .  (3.1) 

The  generator  of  this  bivariate  process  is  easy  to  compute 
since  Vn(t)  subordinates  Xn(t).  Let  f(x,k)  be  a  function 
with  domain  [0,°°)  x  {0,1}  which  is  smooth  as  a  function  of  x 
for  each  k  =  0,1.  The  generator  is  given  by 
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Anf (x,k)=  \ 


'6n[f(x  +  l//n,  k)  -  f(x,k)] 

+  (l-k)nn[f(x  -  l//n,k)  -  f(x,k)]  +  nQf(x,k) 

for  x  l//n 

f 5  n [f (l//n,  k)  -  f(0,k)]  +  nQf(0,k) 

for  x  =  0  . 


(3.2) 


where  Q  is  the  generator  of  the  M/M/1/1  voice  process  and  is 
given  by  (2.1)  . 

Next  expand  in  a  Taylor  series  and  collect  terms  to 
rewrite  (3.2)  as 


Anf(x,k)  = 


nQ f ( x , k )  +  /n  fx(x,k)  (6  -  (l-k)n) 

+  |  fxx(x,k)  («  +  (l-k)n)  +  o(l) 


for  x  >  l//~n 


nQf(0,k)  +  /n  fx(0,k)5  +  i-  fxx(x,k)6  +  o(l) 


for  x  =  0 


(3.3) 


where  f  (x,k) 
x 


f(x,k)  and  fxx(X'k>  =  —  f(x,k) 


We  wish  to  let  n  and  focus  attention  on  the  data 

queue  process  alone.  To  derive  a  diffusion  limit  for  the  data 
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queue  process  alone,  we  introduce  a  sequence  of  functions 

<f  (x,k)  >°°  .  for  which  f  (x,k)  -*■  f(x),  where  f  is  a  smooth 
n  n=l  n 

function  satisfying  f 1 (0)  =0.  We  wish  to  study  the  limiting 

behavior  of  A  f  (x,k).  Let 
n  n 

fn(x,k)  =  f ( x)  +  -i-  g(x,k)  +  ^  h(x,k) 

/n 

where  g  and  h  are  smooth.  Clearly  fn(x,k)+  f(x).  Sub¬ 
stitution  in  (3.3)  gives 

nQf(x)  +  /n[Qg(x,k)  +  f '  (x)  (S  -  (1-k)  n )  ] 

+  [Qh  (x,  k)  +  gx(x,k)  (6 -(l-k)n)  +  |  £' "  (x)  (6  +  (1-k)  n )  ] 

x  >  1//H  (3.4) 

nQf  (0)  +  /JT  [Qg  ( 0  ,k)  +  f '  (0)5  ] 

+  [Qh (0 ,k)  +  g  (0,k)6  +  i  f"(0)(5)]  x  =  0  . 

X 

We  first  note  that  Qf(x)  =  0  since  Q  operates  on  the 
voice  or  k  component  only,  not  on  the  data.  This  eliminates  the 
first  terms  nQf(x)  and  nQf(0).  We  next  examine  the 
Qg(x,k)  +  f '  (x)  (5  -  (l-k)n)  term.  This  can  be  rewritten  as 
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Qg  (x,  k)  +  f'  (x)  [6 


-  (l-k)n  -  ( (<5  -n )  +  <5  pv)/(1  +  pv) 


+  f'(x)((6-n)  +<Spv)/(i  +  pv) 


=  Qg (x, k)  +  f ' (x) 


pv 

-  rrV  +  k  = 

L  v 


f ' (x)  n (l-p) 


We  now  select  the  function  g(x,k)  so  that 


Qg  (x,  k)  =  -f '  (x)  -  y~+  +  k  n 


The  function  g(x,k)  must  satisfy  the  equations 


X (g (x, 1)  -  g(x, 0) )  =  f’(x) 


-u(g(x,i)-  g(x,0) )  =  -f '  (x) 


(3.5) 


The  above  equations  are  consistent  and  redundant,  thus 


any  g(x,k)  for  which 


g (x, 1 )  -  g (x, 0)  = 


f'  ( x) r 

yd  +  p. 


will  suffice.  We  select 


g (x, 0 )  =  \  f'  (x) 


g(x,l)  =  (i  + 


Ud  +  Pv) 


-)  f’(x) 


(3.6) 


9(x'k)  +  <7  +  OT  f0-t>  f'<*> 
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This  choice  of  g  gives 


Vn(,'kl 


-v'  n  n  (1-p)  fMx)  +  [Q  h  { x ,  k )  +  I  ^  1  f "  (x)  (6  -  (1-k)  n) 


+  i  f"  (x)  (6  +  (l-k)n)  1 , 


nk_ 

V' 

x  >  l//n 


(3.7) 


n (pv-k (l+pv) 


/n  If*  (0)6  +  f 1  (0 )  - 


Qh(0'kl  +  (  I  +  vTTfpp" )  f ” 


(0)6  +  j  f" (0)6 


x  =  0 


Equation  (3.7)  can  be  rewritten  recalling  f’(0)  =0  as 


A  f  (x,k)  = 
n  n 


-/n  nd*  p)  f  1  (x)  + 


Qh  (x,k)  +  nf"(x)  J  \  +  -(ln^p-T  (pd-(l-k)) 


+  I  (Pd  +  (1-W)  j 


+  o  (1)  ,  x  >  l//n 


(3.8) 


Qh(0,k)  +  nf"  (0) 


l  ,  nk  .1 

-  +  — m - r  +  o  p, 


pd  2  ^  u  (l+pv)  +  2  pd 


,  x  =  0 


or 


A  f  (x,k) 
n  n 


-/n  nd-p)  f '  (*)  + 


Qh  (x,k)  +  f"  (x)  n*  Pd  I  1  +  y-(l+p  <  11  +  0(1] 

'  \ 


^7) 


x  >  l//n 


(3.9) 


Qh  (0 ,  k)  +  f"  (0)  np. 


(  1  +  T7TTT77  \  )  +  o(1)  ' 


x  =  0 


We  rewrite 


Qh (x,k)  +  f"(x)npd 


as 


4>h  (x,k)  +  f  '  (x) 

and  choose  h(x,k) 


r|2pd  (v. 

u  (L+pv>  V 

so  that 


f"(x)npd 


npy 


Qh  <x,k) 


-f" (x) 


The  function  h(x,k)  must  satisfy 


A  (h  (x,  1)  -  h  (x,0) ) 


f  "  (x) 


n  pdpv^ 
Md+Pv): 


-U(h(x,l)  -  h(x,0) ) 


f"  (x) 


n  p. 


y  (1+pv) 


(3.1*0) 


These  equations  are  consistent  and  redundant  so  a  manifold 
of  solutions  are  possible.  Any  one  will  suffice.  This  choice  of 
h(x,k)  allows  one  to  rewrite  (3.9)  as 


23 


-  t'ii  n(l-p)  f '  (x)  +  f  "  (x) 


np. 


HP, 


vy(i+pv) 


1  + 


+  0(1) 


A  f  (x,k) 
n  n 


x  >  l//ri 


f"  (0) 


o  (  nPv 

dVJ(1+pv) 


1 


+  0(1) 


(3.11) 


We  now  let  n  In  order  for  (3.11)  to  converge  to  a 

sensible  limit,  we  must  invoke  the  heavy  traffic  approximation, 
that  p  ■*  1 .  Specifically  we  let  p  =  p  =  1  -  (0//n)  for  some 
0  >  0.  Equation  (3.11)  becomes  in  the  limit 


A  f  (x,k)  - > 

n  n  r 


Af (x)  =  J 


-n0  f ' (x)  +  f " (x) np 


np. 


d  Vm(1  +  Pv)2 


+  1  » 


f"  (0) np. 


V  (0)  =  0 


np. 


,y(1+Pv) 


T  +  1  ' 


x  >  0 


x  =  0  (3.12) 


This  is  the  generator  of  a  Wiener  process  with  reflecting 
barrier  at  0,  drift  of  -n^,  and 


np. 


T  =  npdl 


^d+Pv) 


+  i 


The  stationary  distribution  for  such  a  process  is  an  exponential 
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distribution  with  parameter 


+ 


r 


j 

thus  the  mean  queue  length  is  given  by 


HP, 


^d+Pv> 


It  is  informative  to  try  to  apply  this  heavy  traffic 
result  to  a  case  in  which  p  is  near  but  less  than  1.  One  might 
replace  0  by  /n  (1-p)  .  Since  we  have  sealed  by  /n,  the 
stationary  distribution  of  N(t),  the  unsealed  queue  length 
should  be  approximately  exponential  with  parameter 

(1-P) 

/  np 

l  v 

d  V<i+Pv>* 


The  mean  queue  length  becomes 


HP, 


1_p  Vy  d+p„) 2 


+  l 


This  result  is  in  exact  agreement  with  Fischer  (1977)  and 

is  therefore  exact.  The  fluid  flow  approximation  treats  n/y  as 

2 

being  large,  hence  [np^/y  (1+p v)  ]  is  assumed  to  dominate  1.  This 

2 

gives  (n  PdPv) / [y (i+Pv)  (1-p)]  as  the  fluid  flow  mean  queue  length 
found  in  Gaver  and  Lehoczky  (1979)  and  in  Section  2  of  this  paper. 
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The  diffusion  approximation  thus  found  represents  a  refinement  of 
the  fluid  flow  diffusion  approximation.  Even  if  n/y  is  large, 
Pv  may  be  small,  so  the  1  may  be  important. 

Once  the  Wiener  process  infinitesimal  drift  and  variance 
have  been  found,  one  can  also  use  the  dyanmics  of  the  Wiener 
process  to  model  the  dynamics  of  the  data  queue  process.  Busy 
period  distributions,  areas,  and  transient  behavior  in  general 


can  be  determined. 
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